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ABSTRACT 

The dynamics of linear perturbations is studied in magnetized plasma shear 
flows with a constant shearing rate and with gravity-induced stratification. The 
general set of linearized equations is derived and the two-dimensional case is 
considered in detail. The Boussinesq approximation is used in order to exam- 
ine relatively small-scale perturbations of low-frequency modes: Cravito-Alfven 
waves (GAW) and Entropy Mode (EM) perturbations. It is shown that for flows 
with arbitrary shearing rate there exists a finite time interval of non-adiabatic 
evolution of the perturbations. The non-adiabatic behavior manifests itself in a 
twofold way, viz. by the over-reflection of the GAWs and by the generation of 
GAWs from EM perturbations. It is shown that these phenomena act as efficient 
transformers of the equilibrium flow energy into the energy of the perturbations 
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for moderate and high shearing rate solar plasma flows. Efficient generation of 
GAW by EM takes place for shearing rates about an order of magnitude smaller 
than necessary for development of a shear instability. The latter fact could have 
important consequences for the problem of angular momentum redistribution 
within the Sun and solar-type stars. 

Subject headings: MHD — Sun: magnetic fields — Sun: rotation — waves — 
stars: rotation 



Introduction 



It is well-known that the most advanced stellar models, taking into account rotation- 
induced hydrodynamic processes of meridional cir culation and shear mixin g, are quite ade- 
quately describing the structure of massive stars (IMaeder fc Meynet II2OOOI ). For solar-type 
stars with relatively wide convective regions, however, these models predict large rotation 
gradients and are in a notab le contradiction with recent observational helioseismology data 
( ICharbonnel &: Talon II2005I ). This is probably due to the insufficiently high efficiency of the 
common hydrodynamic instabilities that are characteristic to these models: these instabil- 
ities can redistribute the angular momentum between the different layers of the radiative 
region, but they are not strong enough to ensure the onset of a uniform rotation regime. 

While the radiative interiors of solar-type stars rotate quite uniformly, their convective 
zones rotate differentially. Mounting observational evidence indicates that the radiative 
interior of the S un, for instance, rotates almost as a solid body, with a quasi- flat seismic 
rotation proffie (IChaplin et al. Ill999l : ICouvidat et al. II2003I ). But the solar convective zone 
exhibits a strong shearing in the latitu dinal direction, with its e quatorial layers rotating 
about 1/3 faster than the polar regions (IKim fc MacGregor 1120031 ). The transition between 
the differential and the uniform rotation regimes takes place in a relatively thin layer (with 
thickness < 0.05-Rq), called the tachodine, where t he radial shearing of the rotation is 
particularly significant (IMiesch II2005I : iPetrovay 1120031 ). It is often argued that some kind of 
shear instability must be the s ource of gravity waves at the base of the solar convection zone 
( iKumar. Talon &: Zahn Ill999l ). either in the immediate vicinity or just inside the tachocline. 
This instability can hardly be of a purely hydrodynamical nature, because from various 
numerical studies of the buoyant rise of magnetic flux from the bottom of the convection 
zone to the surface of the sun it follows th at the tachocl ine is strongly magnetized with its 
poloidal magnetic field about 10^ - 10^ G jMiesch hopj ). 



From this perspective it seems quite plausible to surmise that on the time-scale of the 
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order of the solar age or less not only purely hydro dynamic internal gravity waves (IGWs) 
participate in the process of the angular momentum redistribution within the solar interior, 
but also , or rather, Gravito-MHD waves. This idea is not new. I n the literature it has been 
argued (jZahn. Talon &: Matias 1 119971 : iKumar &: Quataert 1 119971 ) that Gravito-MHD waves 
could be one of the most promising candidates for the onset of the redistribution processes. 
It was argued that a self-consistent model should comprise a large-scale magnetic field in the 
Sun's interior, and an accurate consideration of the Coriolis effects in the convection zone 



and in the tachocline (ITalon fc Charbonnel 1 1 19981 . l200l 12004 l2005l : ICharbonnel fc Talon 
20051 ) . It was also argued that turbulent stresses in the convection zone, via the agency of 



Coriolis effects, induce a meridional circulation, causing the gas from the convection zone 
to burrow downwards and generate both horizontal and vertical velocity shear characteriz- 
ing the tachocline. However, the interior magnetic field confines the burrowing and, hence, 
smooths and diminishes the shear, a s dema nded by the observed structure of the tachocline. 
According to I Charbonnel fc Talon I (120051 ) . a decisive test of this qualitative scenario, af- 
ter numerical refinements, would be its quantitative consistency with the observed interior 
angular velocity structure. 

The aim of the present paper is to present the results of a detailed study of the hnear 
dynamics of perturbations in gravitationally stratified magnetized shear flows. Bearing in 
mind the highly probable importance of the Gravito-MHD waves for the structure of the Sun 
and solar-type stars, we decided to pose and study the following three interrelated issues: 
1) how the presence of the nonuniform velocity field affects the propagation of the waves 
through the stellar plasma; 2) what kind of energy exchange processes between the different 
collective modes and between the modes and the ambient flow may happen; and 3) what 
other astrophysical consequences these processes could have. 

Originally, we develop the three-dimensional (3D) model, allowing the ambient flow to 
have velocity shearing in both transversal (prior to the gravity field) directions. Further on, 
the study is focused by a number of simplifying assumptions, enabling the solution of the 
equations. First of all, we consider only two spatial dimensions (2D) bearing in mind that it 
is quite straightforward to extend the analysis to the fully three-dimensional case. Second, 
we assume a constant steady-state (or s tationary) shear flow, an assump tion that allows us 
to use the shearing sheet approximation (iGoldreich &: Lynden-Bell Ill965l ). And third, we use 
the Boussinesq approximation, and hence we study the dynamics of relatively small-seal^ 
perturbations of the low-frequency modes, viz. the Gravito-Alfven waves (GAW) and the 
Entropy Mode (EM) perturbations. 
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Our study leads to the conclusion that for arbitrary shear in the flow, there exists a 
finite interval of time in which the evolution of the modes is highly non-adiabatic. Outside 
this interval, the evolution is adiabatic and can be accurately described by a WKB approx- 
imation. Therefore, we have an asymptotic problem, quite common for various quantum- 
mechanica l applications, where one needs to obtain connection formulae in turning point 



problems ( iLandau fc Lifshitz 1 119771 ) . The relevant analysis is done in this paper and it is 
found that the non-adiabatic behavior of the considered modes manifests itself in the form 
of two phenomena, viz. the over-reflection of the GAW and the generation of the GAW by 
the EM perturbations. We show that for flows with moderate and high shearing rates both 
these processes are very effective converters of the equilibrium flow energy into the energy 
of the waves. 

In the context of above-discussed astrophysical problems our results imply that efficient 
generation of the GAW by the EM perturbations takes place for shearing rates of about an 
order of magnitude smaller then necessary for development of shear instability. 

The remainder of the paper is organized in the following way: the general mathematical 
and heuristic formalism is presented in the Sec. 2. The over-reflection of GAW is studied in 
the Sec. 3; and the second non-adiabatic process, viz. the generation of the GAW by the EM 
perturbations, is studied in the Sec. 4. Finally, in Sec. 5 we conclude with a summary and 
a brief discussion of our results in the context of their possible importance for the angular 
momentum redistribution within solar-type stars and the solution of the problem of the 
uniform rotation of their radiative zones. 



2. Linear perturbations of sheared plasma flows 

For the sake of generality, we derive the basic set of equations for the case of 3D per- 
turbations in a compressible, gravitationally stratifled MHD fluid. This formalism will be 
presented in the next (2.1) subsection. Later on (subsection 2.2), we will restrict the con- 
sideration to the somewhat simpler case of 2D perturbations in an incompressible medium. 
Besides, the Boussinesq approximation will be used and the stationary shear flow will be 
considered. 



2.1. General theory 



In our model, the geometry of the considered problem is simplifled in the following way: 
the equilibrium flow U is supposed to be plane-parallel, to be directed along the x-axis, and 
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to have both a horizontal (Ay) and a vertical (A^) shear: 

V=[Ayy + A,z,0,0]. (1) 

The uniform gravity g is assumed to be directed along the negative direction of the 
z-axis: 

g = [0, 0, -g] . (2) 

We consider a simplified model and assume that the equilibrium magnetic field Bq is 
toroidal, parallel to U and that it is possessing the gravity-induced vertical stratification: 

Bo= [B{z), 0, 0]. (3) 

The set of equations of one- fluid ideal magnetohydrodynamics (MHD), governing the 
physics of the flow is: 

Ap + p(V-V) = 0, (4) 
pDtV ^ -v(^P + + ^{B ■V)B + pg, (5) 

AB= (B- V)V-B(V- V), (6) 

DtS = 0, (7) 

while 

V • B = 0. (8) 

In (4-7) the total (convective) time derivative operator is denoted by: Dt = dt + {V -V). 

Instantaneous values of all physical variables are decomposed into their regular (mean) 
and perturbational components: 

P = Po + Q, (9) 

P = Po+P: (10) 

S = So + s, (11c) 

V = U + u, (12a) 

B = Bo + b. (126) 

It is straightforward to see that the horizontal background flow does not affect the 
vertical hydromagnetic equilibrium of the flowing MHD fluid, governed by the equation: 

d,(^Po + ^^=-po9. (13) 
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In general the variation of the density dp is related with variations of the pressure dP 
and entropy ds via the relation: 



where: 



dp = pds + {1/Cf)dp, 

p = {dp/ds)p, 
Cl ^ (dP/dp),. 



(14) 

(15) 
(16) 



Linearized set of equations governing the evolution of perturbations within this flow can 
be written in the following way [Vt = dt + {AyU + Azz)dx]: 



'DtQ + {dzPo)v^ + po(V • u) = 0, 



Alia 



—Auy — Auz dxP + 

Po 



VtUy = - — dyP - ( ) {dj)y - 

Po v47rpo/ 



47rpo 



(17) 
(18) 

(19) 



— ( 7^ ]{dxhz 



h 3 

Ox Q, 

Po 



47rpo 

Vts + {dzs)vz = 0, 
Vthy = BodxUy, 
Vfbz = BodxUz, 

dxbx + dyby + d^b^ = 0. 



(20) 

(21) 
(22) 
(23) 

(24) 



Due to the stratification some coefficients on the right hand sides of Eqs. (17-23) depend 
on z. But when studying the dynamics of small-scale perturbations (with <^ Zq) one 
can consider the mean components as constant. In this case the set of first-order, partial 
differential equations (17-24) can be reduced to the set of ordinary differential equations 
(ODEs) with time-dependent coefficients if we look for solutions in the following form: 



with 



F{r,t) = F[k(t),i]e'[('^")-^('^'*)1, 
^(k,t)= (v- [k{t)dt], 



(25) 
(26) 
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where k(t) is a time-dependent wavenumber vector, determined by the following s et of equa- 
tions (iLagnado et al. 1 1 1984 ICraik fc Criminale I Il986l : iMahajan fc Rogava l ll999h : 



k«+5^.k = 0, 



(27) 



with being the transposed shear matrix: 



Ux.x Ux y Ux 



a. 



y,x 



^y,y ^y,z i ) 



(28) 



z,y 



which, in our case has only two nonzero components: Ay and A-,: 



Ay Az 







(29) 



Due to Eq. (27), transversal components of the \s.{t) acquire linear time-dependence: 

ky{t) = ky{0) - Aytkx, (30a) 

hit) = kM - Aztkx. (306) 

Therefore, we see that the wave-number vector k()f:) of the Spatial Fourier Harmonics 
(SFH) varies in time, i.e. in the linear approximation there is a "drift" of the SFH in the 
k-space. The physical reason of this "drift" is related to the fact that in the sheared flow the 
perturbations cannot have the form of a simple plane wave due to the effect of the shearing 
background on the wave crests. 

Applying the ansatz (25) to Eqs. (17-23), we reduce the system to the following set of 
first-order ODEs [hereafter we are using the following notation /*^") = d^f]: 



o?(i) = ievz + Vx + ICy{T)vy + /C^(r)f, 



V 



(1) 



-TZvy — RVz 



V, 
V 



(1) 

y 

(1) . 



-CJ^lCy{T)V + by- ICy{T)bx 

-a^lCz{T)V + K-ICz{T)bx- 



5(1) 
y 



yi 



(31) 
(32) 
(33) 

(34) 

(35) 
(36) 
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6« = -v.. (37) 
while the Maxwell equation (24) reduces to the following algebraic relation: 

6,. + /Cj,(r)6j, + /C,(r)6, = 0. (38) 

In these equations, the following dimensionless notation is used for the constants: 
Ry = Ay/KCA, R, = AJk.CA, }Cy{0) = ky{0)/k,, /C,(0) = k,{0)/k,, e = {Kz^y\ = 
{CsICaY, N"^ = {Nbv /CAkxY ■ The dimensionless variables appearing in the above set 
of equations are defined as: Vj = Uj/CA, bj = ibj/Bo, d = ig/po, V = ip/poCg, e = 
-kJ/{d,so), T = kxCAt, )Cy{T) = /C,(0) - 7^r, /C,(r) = /C,(0) - Rt. 

Notice that this system is not closed, because we have only eight equations for nine 
variables (v, b, d, e, V). The closure of the set of equations is ensured by the thermodynamic 
relation that follows from Eq. (14) [a = g/k^C"^]: 

d = V -i{N'^/a)e. (39) 

This system of equations describes the temporal evolution of Gravito-MHD waves mod- 
ified by the presence of a velocity shear in the two planes transversal to the flow direction. 
The full analysis of this set of equations is beyond the scope of the present paper. Instead, 
in the next sections, we focus our study on the relatively simple 2D and incompressible 
case. We will see that even in this simplified case the presence of the velocity shear brings 
a considerable novelty in the dynamics of perturbations. 

2.2. The incompressible, 2D limit 

In order to simplify the mathematical aspects of the problem, we make the following 
assumptions: 

• We consider the 2D case, restricting the problem to the xOz plan and implying: 

}Cy{0) =Ry = 0. (40) 

Hereafter we write R = Rz{A = Az) In this case, Eq. (38) enables us to express the lon- 
gitudinal (x-) component of the magnetic field perturbation in terms of its transversal 
(2;-) component: 

bx = ~ICz{r)bz. (41) 



^The Alfven speed — B^/Attpq and the Brunt- Vaisala frequency N^y — —gdzPo/po are defined in the 
usual way. 
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We drop the S^^^ term in th e continuity equatiori, i.e. we adopt the concept of so-called 
dynamic incompressibility (iLandau fc Lifshitz Ill959l ). This implies that the velocity 
of the flow is assumed to be small enough in comparison with the speed of sound. 
As a consequence, the density perturbations evoked by the pressure perturbations are 
negligible and instead of the Eq. (39) we will have: 



d 



-i{N'^/a)e. 



(42) 



However, since the medium is thermally conducting its density should vary also due to 
the temperature variation, and this effect should be (and is) taken into account. 



In hydrodynamics, the second assumption is known as the Boussinesq approximation: in 
the conservation equations the terms proportional to the density perturbation are retained if, 
and onl y if, they produce a buoyancy force § . The Boussinesq approximation automatically 
imphes dhandau fc Lifshitz 1^9591 : Icill Ill982h that the vertical length scale of perturbations, 
^z ~ is much smaller then the characteristic stratification length scale zq. This condition 
was already used for the derivation of Eqs. (31-37). In terms of our dimensionless notation, 
this condition transcribes into: 

£ < 1. (43) 

When this condition is satisfied, the above approach enables us to study the dynamics of 
small-scale, low-frequency perturbations. 

These assumptions make the basic set of equations much simpler. Taking into account 
(41) and (42), we get [with 10{t) = 1 + JCKr)]: 



Vx + KLz{t)vz = 0, 



(44) 

(45) 
(46) 
(47) 
(48) 



^The last term in Eq. (34). 



"^ for detailed analysi s of th e Boussinesq approximation in stratified shear flows see iDidebulidze I (|l997l ) 
and bidebulidze et. all |2004l ). 
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From the last two equations it is apparent that there is an algebraic relation between 
the entropy perturbations and the magnetic field perturbations, viz. 



e + 62 = X = constir). 



(49) 



Another, less obvious, relation between the perturbed quantities may be found if one 
takes time derivative of (44) and takes into account (45) and (46). The result can be used 
to express the pressure perturbation as a function of the other variables, i.e. 



-2R W^P 
-Vy + (Jby - , e. 



(50) 



Having three algebraic relations, Eqs. (44), (49), and (50), for five variables, we can 
reduce the system to one second-order ODE. For instance, for the variable &(r), defined as 
6(r) = }C{T)bz, we can derive the following explicit, second-order inhomogeneous ODE: 



dr2 



+ 



1 + 



iV2 



/C2(r) IC' 



-X. 



(51) 



This equation descr i bes th e dynamics of the so-called continuous eigenspectrum (see, e.g. 
Balmforth fc Morrison I (120021 ) and the references therein). In order to take into consideration 
the perturbations belonging to the discrete eigenspectrum one should specify some boundary 
conditions with respect to the vertical spatial variable z. In the absence of an external 
magnetic field, it is known (IDrazin fc Reid Ill98ll ) that all the discrete eigenmodes are stable 
if the Richardson number, defined as Ri = N'^y/A^ = N'^/R'^, is greater then 1/4. 

Equation (51) describes two important physical phenomena related to the dynamics 
of the continuous eigenspectrum, which will be considered in detail in the two subsequent 
sections, viz. (a) the (over)-refiection of the GAW and (b) the conversion of the EM into tghe 
GAW. It will be shown that the EM perturbations can be a very efficient source of GAW, 
even for relatively small shearing rates {Ri > 1/4). 



Over-reflection of Gravito-Alfven waves 



In order to describe the phenomenon of the over-reflection of the GAW in its pure form, 
let us study Eq. (51) with T = 0, yielding: 



fb 



Ar2 



R 



2 1 



0. 



(52) 
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The standard meth ods of the analysis of simi lar equations (jOlver Ill974j ) is well-known 
in quantum mechanics (ILandau &: Lifshitz 1 119771 ) . The behavior of the solution is fully 
determined by the shear-modified frequency: 



i?2 



JCHt) JCHt)' 



(53) 



From the mathematical point of view, the (over)-reflection is caused by the cavity of 
00"^ (t) that appears in the vicinity of the point t*, where Kz{t*) = 0. 

In the areas where Kz{t) ^ R, the evolution of 6(r) is adiabatic [dct;(r)/dr ^ ^^^(t)] 
and, therefore, the approximate solution has the form 



6(r) ^ A+6+(r) + A_6_(r) 



where 



(54) 



(55) 



are WKB solutions having positive and negative phase velocities along the x-axis, respec- 
tively, while the A± can be treated as the intensities of the corresponding types of pertur- 
bations. 

Let us assume that and A^ are the WKB am plitudes far o n the left and right hand 
side from the point r*, respectively. It is well-known (jOlver Ill974l ) that the conservation of 
the Wronskian leads to the following invariant: 



\A 



L |2 



I A 



L |2 



\A 



R\2 



\A 



R\2 



(56) 



which physically corresponds to the conservation of the wave action (see, e.g.. lGogoberidze et al. 
(120041 ) and references therein). If initially A^ = 0, the refraction and transmission coeffi- 
cients may be defined in the usual way: 



Re 



|A^P 



\A 



L 12 ' 



D 



(57) 



while from Eq. (56) we obviously have: 



l + Re = D. 



(58) 



The dependence of the transmission coefficient D{R) on R, obtained through the nu- 
merical solution of Eq. (52) for different values of the Brunt- Vaisala frequency, viz. = 
(solid line), = 1 (dashed line) and N = 2 (dash-dotted line), is presented on Fig. 1. The 
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Fig. 1. — Dependence of the transmission coefficient D{R) on tfie normalized sfiear parameter 
R for different values of the Brunt- Vaisala frequency: = (solid line), = 1 (dashed 
line) and N — 2 (dash-dotted hne) . 
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initial conditions for the numerical solution are chosen from the WKB solutions (55) far on 
the left-hand side of the resonant area (-^'^(O) ^ 1, -R)- The analysis leads to the conclusion 
that the amplification of the energy density of the perturbations is always finite, but it can 
become arbitrarily large with a proper increase of the shearing rate. It can be seen that the 
velocity shear can cause a strong amplification of the GAW, even for moderate values of the 
normalized shearing rates. 

According to the numerical study of Eq. (52), the amplitude of the reflected wave exceeds 
the amplitude of the incident wave if: 



> 2. 



(59) 



This phenomenon, originally discovered by Miles for acoustic waves (iMiled 119571 ). is called 
over-reflection. 

Bearing in mind that the energy of the mode is proportional to the square of its ampli- 
tude, we can easily surmise that in the course of its non-adiabatic evolution the perturbation 
energy is increasing. According to Eqs. (56) and (58) the ratio of the total 'post-reflection' 
energy to the initial one equals Re + D > 1. The energy increases at the expense of the mean 
(shear) flow energy. Consequently, over-reflection represents a quite efficient mechanism for 
transferring the mean flow energy to perturbations. 

In the limit — > 0, the reflection coefficiei it becomes exponentially small with respect 
to the parameter —l/R, i.e.. Re ~ exp(— 1/i?) (iGogoberidze et al. 1 120041 ). 



4. The generation of Gravito-Alfven waves 

Let us return to the analysis of the full, inhomogeneous equation (51). When the 
condition of the adiabatic evolution da; (r)/dr -C uj'^(t) holds, then the ap proximate solution 
of Eq. (51) has the following form (IChagelishvili. Rogava &: Segal Ill994l ): 

6(r) ^ A+6+(r) + A_6_(r) + 6s(r), (60) 

where the non-periodic solution 6£;(r) is given by: 

Ar2 / i?2 x-i 



and can be readily identified as the shear-modified Entropy Mode (EM) perturbation. 

It turns out that another physically important phenomenon that takes place in the 
system under consideration for moderate and high shearing rates, is the generation of the 
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Fig. 2. — Temporal evolution of 6(t) forN — l and R = 0.7. 
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GAW by the EM perturbations. An illustrative example is given by the Fig. 2, where the 
numerical solution of Eq. (51) for the case when = 1 and R = 0.7 for the variable 6(r) is 
shown. The initial conditions were chosen in such a way that at the beginning (at Ar^(O) = 35) 
there exists only the EM perturbation. Namely, the initial conditions were found by means of 
Eqs. (60) and (61) with A± = and X = 1). The figure clearly demonstrates the generation 
of the GAW while the EM passes through the interval of non-adiabatic evolution, situated 
in the vicinity of the point /C^(r*) = 0. 

For the quantitative description of the wave generation process, analogously to the 
reflection and transmission coefficients defined in the previous section, we define the wave 
generation coefficient G± in the following way: 



G± = 





2 


|X 


2JV2 



With this definition, G± is proportional to the ratio of the generated wave energy to 
the energy of incident entropy mode perturbation. From the symmetry of the problem, it is 
obvious that G+ = G- = G. 

In Figure 3, the dependence of the generation coefficient G{R) on the normalized shear 
parameter R for different values of the Brunt- Vaisala frequency: A^ = 1 (solid line), N = 2 
(dashed line) and N = 3 (dash-dotted line) is presented. 



5. Discussion 

In order to explain the near uniformity of the rotation profile in the radiative region of 
the Sun and similar regions in solar-type stars usually two different mechanisms are being 
proposed. The first o ne imphes the presence of tr avelling IGWs generated at the base of 



the convective region (jCharbonnel fc Talon 1 12005| ). while the second one is based on the 



presumed (and for the Sun observational ly justified) existence of the magnetic field in the 
radiative zone jGough fc Mclntvre lEJ l 



It seems reasonable to assume that the above-mentioned two mechanisms are not nec- 
essarily alternatives, but could be envisaged as complementary ones. In particular, we may 
argue that if magnetic fields are indeed present inside the radiative zone and the Gravito- 
MHD waves do pass through the zone, it leads to intense energy exchange processes between 
the waves and the stellar plasma. This circumstance, in its turn, might lead to the more 
efficient redistribution of the angular momentum within these stars, establishing quasi-fiat 
rotation profiles. The presence of the magnetic field implies that it is more reasonable to 
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Fig. 3. — Dependence of the generation coefficient G{R) on the normahzed shear parameter 
R for different values of the Brunt- Vaisala frequency: N = 1 (solid line), N = 2 (dashed 
line) and N — 3 (dash-dotted fine) . 
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25 




Fig. 4. — Dependence of the generation coefficient G{Ri) on Ricfiardson number Ri for 
AT = 4. 
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speak not about purely hydrodynamical IGWs, but about Gravito-MHD waves in general, 
and Gravito-Alfven waves in particular. 

The evolution of the latter mode, present only in the magnetized case and absent in a 
purely hydrodynamic medium, is studied in the present paper. The theoretical background 
for studying of the evolution of the full set of Gravito-MHD waves is also developed. 



From the general theory of IGWs, we know ( iDrazin &: Reid Ill98ll ) that all the discrete 
eigenmodes are stable under the condition that Ri > 1/4. Our study shows that the genera- 
tion of the GAW by the EM perturbations could provide an efficient mechanism for creating 
the required flux of the GAW, even for relatively low shearing rates, when the discrete eigen- 
modes are stable. Note, in this context, that for the case presented in Fig. 2, for instance, 
the Richardson number was about Ri ^ 1.02. 

In the magnetized case, we studied the dependence of the generation coefficient G{Ri) 
on the Richardson number Ri for the relatively weak external magnetic field (A^ = 4). The 
results of this study are presented in Fig. 4. The results clearly indicate that even though the 
generation coefficient falls rapidly with increasing Ri, the efficient generation of the GAW 
by the EM perturbations still takes place for Ri < 1.25. 

Obviously, our study has a simplified nature and it only indicates at the possibly im- 
portant role of non-modal phenomena related with GAW in the redistribution of the angular 
momentum and the flattening of the rotation curves in solar-type stars. For making this 
claim more convincing, we have to consider the compressible, three-dimensional analogue of 
this problem and we have to verify whether the full set of Gravito-MHD waves is as efficient 
(or even more efficient) in exchanging energy with the equilibrium flow. We know from previ- 
ous studies that in the absence of the gravity-induced stratification among the shear-modified 
MHD waves, the fast mag netosonic waves are the most efficient ones in exchanging energy 
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20001 ). It is of great interest to check whether the fast Gravito-magnetosonic waves have the 
same quality and to determine what their contribution in the angular momentum redistri- 
bution could be. The study of the linear dynamics of all Gravito-MHD waves is currently 
initiated and the results will be published in a subsequent paper. 

Finally, one has to remember that the non-modal approach that has been applied in 
the present paper, provides no information about the spatial aspects of the shear-induced 
processes, because the study of the Spatial Fourier Harmonics is conflned to the phase space 
of the wave number vectors k(t). In order to have a clear idea about the spatial appearance of 
these processes, one has to study them numerically (similarly to the uniform, non-stratifled 



flow study made by lBodo et al. I (120011 )) and to check how the non- modal phenomena couple 
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with the traditional rotational instabilities. 
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